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Abstract 

We discuss the effects of multi-phonon excitations on heavy-ion fusion reactions at 
energies near and below the Coulomb barrier, focusing especially on the role of anhar- 
monicities. We carry out a systematic study of those effects on the excitation function 
of the fusion cross section and on the fusion barrier distribution, by using the vibrational 
limit of the interacting boson model. We also analyze the recently measured high-precision 
data of the 16 O + 148 Sm fusion reaction with this model and discuss the anharmonic prop- 
erties of the quadrupole as well as the octupole vibrations in 148 Sm. Negative and positive 
static quadrupole moments are deduced for the first 2 + and 3~ states in 148 Sm, respec- 
tively. It is shown that the fusion barrier distribution strongly depends on the sign of 
the quadrupole moments, suggesting that subbarrier fusion reactions offer an alternative 
method to extract the static quadrupole moments of phonon states in spherical nuclei. 

PACS number(s): 25.70. Jj, 21.60.Fw, 21.10.Ky, 27.60. +j 
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I. INTRODUCTION 



Large enhancement of heavy-ion fusion cross sections at subbarrier energies relative 
to predictions of a one- dimensional potential model has been well recognized [I]. It 
is by now well established that these enhancements can be attributed to couplings of 
the relative motion between the colliding nuclei to their intrinsic excitations as well as 
transfer reactions || . When the excitation energy of an intrinsic motion is neglected, the 
corresponding channel-coupling gives rise to a distribution of potential barriers [§, f|. || . 
Based on this idea, a method was proposed to extract barrier distributions directly from 
fusion excitation functions by taking the second derivative of the product of the fusion 
cross section and the center of mass energy Ea with respect to E || . It was recently shown 
that the concept of the barrier distribution can be extended also to the cases where the 
intrinsic motion carries a finite excitation energy J7|. Coupled-channels calculations have 
shown that the fusion barrier distribution, i.e. d 2 (Ea)/dE 2 , is very sensitive to the details 
of the channel couplings, while the fusion excitation function itself is rather featureless 
||. In order to deduce meaningful experimental barrier distributions and to observe 
their sensitivity to the couplings, the excitation function of fusion cross sections has to be 
measured with high precision at small energy intervals. Thanks to the recent developments 
in experimental techniques |§, such data are now available for several systems, and they 
have clearly demonstrated that the barrier distribution is indeed a sensitive quantity to the 



channel couplings [ 10 1 . They clearly show the effects of couplings to static deformations 



and associated rotational motions |TU| , O , |T2| , |T3| Pfl , vibrational degrees of freedom 



14, 15], [16], |17], 0, and transfer reactions fLO], [14] [H| [T(| [17], in a way much more 



apparent than in the fusion excitation function itself. These high precision data have 
thus enabled a detailed study of the effects of nuclear collective excitations on fusion 
reactions, and have generated a renewed interest in heavy-ion subbarrier fusion reactions. 

Nuclear surface oscillations with various multipolarities are typical examples where 
the barrier distribution analyses of fusion reactions may be applied to study their de- 
tailed properties. For instance, the barrier distribution analysis of the recently measured 
accurate data on the 58 Ni + 60 Ni fusion reaction has shown evidence for couplings to 
double-phonon states in 58 Ni and 60 Ni [IB]. The barrier distributions were shown to be 



quite sensitive to the number of phonons excited during fusion reactions. To the zeroth 
approximation, these surface vibrations can be described in terms of harmonic oscilla- 
tions, which dictate simple relations among the level energies and the electromagnetic 
transitions between them. Some of the characteristic features of harmonic oscillations are 
summarized as follows: (i) all the levels in a phonon multiplet are degenerate and the 
energy spacing between neighboring multiplets is a constant, (ii) the electric transitions 
between neighboring multiplets are linearly related, e.g. B(E2; if — > (1—2)+) = (1/2) 
B(E2; l\ — > 0i ), (iii) the static quadrupole moment is zero for all the phonon states. 
In realistic nuclei, however, there are residual interactions between phonons, which cause 
deviations from the harmonic oscillator limit |19|]. The anharmonicity leads to a level 
splitting within a multiplet and to modifications in the ratios between various electro- 
magnetic transition strengths. Furthermore, the anharmonicities generate a finite value 
of static quadrupole moment in excited states PU . 
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In even-even nuclei near closed shells, there are many examples of two-phonon triplets 

Though the central position of their 



(0 + ,2 + ,4 + ) of quadrupole surface vibrations [pi 



excitation energies is approximately twice the energy of the first 2 + state, they usually 
exhibit appreciable splitting within the triplet. A theoretical analysis of the anharmonici- 
ties for the quadrupole vibrations was first performed by Brink et al. |Tj|, who related the 
excitation energies of three-phonon quintuplet (0 + , 2 + , 3 + , 4 + , 6 + ) to those of the double- 
phonon triplets, and also gave relations between the electric transition strengths from the 
three- to the two-phonon states and those from the two- to the one-phonon states. These 
relations, however, had not been confirmed until recently because of the sparse experi- 
mental data on three-phonon states. The experimental situation has improved rapidly 
and data on multi-phonon states are now available for several nu- 
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m recent years 

clei. As a consequence, study of multi-phonon states, and especially their anharmonic 

It is worth mentioning that 



51 26, 27 



properties, is attracting much interest 
anharmonic effects are not restricted to low-lying vibrations but play significant roles also 
in multi-phonon excitations of giant resonances in heavy- ion collisions |28j, |29|, |30[| . 

In many even-even nuclei near closed shells, the first 3~ state appears at a relatively 
low excitation energy, which competes with the quadrupole state ^FL, |31| . These low-lying 
3 _ states have been frequently interpreted as collective octupole vibrations arising from a 
coherent sum of one-particle one-hole excitations between single particle orbitals differing 
by three units of orbital angular momentum 21, [32|. This picture is supported by the large 
E3 transition probabilities from the first 3~ state to the ground state, and suggests the 
possibility of multi-octupole-phonon excitations. In contrast to the quadrupole vibrations, 
however, so far there are only little experimental evidence for double-octupole-phonon 
states. One reason for this is that E3 transitions from two-phonon states to a single- 
phonon state compete against lower multipolarity El transitions. This makes it difficult 
to unambiguously identify the two-phonon quartet states (0 + , 2 + , 4 + , 6 + ). Consequently 
despite the fact that the first 3~ state of 208 Pb has a large quadrupole moment, which 
indicates a significant anharmonic effect in octupole vibrations |53[ |34|, |35], |36|, a direct 
study of the anharmonic properties in multi-octupole phonon spectra has not been possible 
for a long time. Only in recent years, convincing evidence have been reported for double- 
octupole-phonon states, as well as double-phonon states built from single octupole and 
quadrupole phonon states, in some nuclei, e.g., 208 Pb 



37|, 144 Sm P 
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Sm 
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g, 145 Nd [0, 144 Nd H, and 148 Gd g§ |5 
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we have shown that heavy-ion subbarrier fusion reactions provide a pow- 
erful method to study the anharmonic properties of quadrupole as well as octupole vi- 
brational excitations. We analyzed the recently measured high precision data for the 16 
+ 144 Sm reaction and extracted negative quadrupole moments for both the first 2 + and 
3~ states in 144 Sm. The aim of the present paper is, in addition to giving the full de- 
tails of that analysis, to carry out a systematic study of the effects of anharmonicities on 
subbarrier fusion, and to apply the same analysis to the 16 + 148 Sm fusion reaction to 
study the anharmonic properties of the vibrational excitations in 148 Sm and their effects 
on the fusion reaction. Although the importance of the anharmonic effects on proton 
scattering [ 4"7|] , as well as the reorientation effects in phonon spectra on heavy- ion elastic 
and inelastic scattering ESI HJ has been pointed out, no systematic studies of the effects 
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of anharmonicities on subbarrier fusion and on fusion barrier distributions have been per- 
formed so far. In view of the high precision data on subbarrier fusion reactions, which 
have recently become available, such studies are necessary in discussing the effects of vi- 
brational excitations on subbarrier fusion reactions. The paper is organized as follows. In 
Sec. II, we formulate coupled-channels calculations which explicitly take into account the 
anharmonic properties of vibrational excitations using the vibrational limit of the inter- 
acting boson model |50| . We compare the formalism with that in the harmonic limit and 



point out several important features of the anharmonic vibrational excitations. In Sec. Ill, 
a systematic study of the effects of anharmonicities on fusion reactions is presented. They 
include the effects of unequal spacing of levels, reorientation, and finite boson number. 
In Sec. IV, the formalism is applied to the 16 O + 148 Sm fusion reaction. We extract the 
quadrupole moments of the first 2 + and 3~ states in 148 Sm from the analysis of the high 
quality fusion data available for this system. Finally, a summary of the work is given in 
Sec. V. 

II. EFFECTS OF PHONON EXCITATIONS ON FUSION 

In this section, we describe the basic formalism for treating the vibrational excitation 
in subbarrier fusion reactions. Let us first consider the case where the relative motion 
between the colliding nuclei couples to the quadrupole vibrations in the target nucleus. 
The total Hamiltonian of the system is assumed to be 

H = -^-V 2 + V N (r) + ZpZre2 + H int + V coup (r,0, (1) 
l\i r 

where r is the coordinate of the relative motion between the projectile and the target, and 
fi is the reduced mass. Vn is the bare nuclear potential, Z P and Z T are the atomic numbers 
of the projectile and the target, respectively. H int describes the vibrational excitations in 
the target nucleus, and V coup , the coupling between these excitations (generically denoted 
by £) and the relative motion. In Sub-sections A and B, we discuss the coupling of the 
harmonic and anharmonic vibrators within the linear coupling approximation. Although 
this approximation is too simplistic to apply to realistic systems |5T], [52|, it enables us to 



easily understand the essential features of the effects of anharmonicities. Extensions of 
the model so as to include the couplings to all orders and to the octupole vibrations in 
the target nucleus as well are given in Sub-sections C and D, respectively. 

A. Harmonic limit 

In the geometrical model of Bohr and Mottelson, the radius of the vibrating target 
nucleus is parameterized as 



R{6, 0) = Rr 1 + £ oc^Y*^, 0) , (2) 



where R? is the equivalent sharp surface radius of the target nucleus and a 2/1 are the 
coordinates of the surface vibration. To the lowest order, the surface oscillations are 
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approximated by a harmonic oscillator and the corresponding Hamiltonian is given by 



H int = ha; 4 /t «2 M + 2 j • ( 3 ) 

Here hu is the oscillator quanta and aJ>„ and a 2tli are the phonon creation and annihila- 
tion operators, respectively. The surface coordinates a 2fl and the phonon creation and 
annihilation operators are related by 



®2l 



"o (4 M + H"a2 M ) , (4) 



where is the amplitude of the zero point motion. It is related to the quadrupole 
deformation parameter by cto = /^/vS |H] and can be estimated from the experimental 
transition probability using 



1 4vr B(E2) T 

In the collective model, the coupling Hamiltonian between the relative motion and the 
quadrupole surface oscillations is often assumed to be 

V coup (r, a 2 ) = f{r) ^Y^{r). (6) 

This is the so called linear coupling approximation where only the linear term with respect 
to the surface coordinate is taken into account. The coupling form factor, f(r), consists 
of the nuclear and the Coulomb parts and is given by 

/(r) = _ flr ^ + jW^f. (7) 

The dimension of the resulting coupled-channels calculations is in general too large 
for practical purposes. Here, we introduce two approximations which considerably reduce 
the dimension of the coupled-channels calculations with negligible effects on its accuracy. 
The first one is the no-Coriolis approximation where one transforms the whole system 
to the rotating frame where the z axis is along the direction of the relative motion r 
at every instant |53| , |54| . This is achieved by replacing the angular momentum of the 



relative motion in each channel by the total angular momentum J. As the operator 
for a rotational coordinate transformation in the whole space then commutes with the 
centrifugal operator for the relative motion, one can make this transformation to the 
rotating frame without introducing any complications |54|. Using Y 2fJi {j = 0) = JB/^kS^q, 
the coupling Hamiltonian (^) in the rotating frame becomes 

V coup (r, a 2 ) = J^-f(r)a 20 = ~j=f( r ) (4o + a 2o) ■ (8) 



Since the projection of the total angular momentum J onto the z-axis of the rotating 
frame, i.e. the K quantum number, is conserved in this approximation, the dimension of 
the coupled-channels equations is drastically reduced. 
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A further reduction of the coupled-channels equations can be achieved by introducing 
the n-phonon channels [[53], |55| . For the quadrupole surface vibrations, the two phonon 
state has three levels (0 + ,2 + ,4 + ). In the harmonic limit, this two-phonon triplet is de- 
generate in the excitation energy. One can then replace the couplings to all the members 
of the two-phonon triplet by the coupling to a single state given by 

|2 >= < 2020|/0 > |/0 >= ^(4 ) 2 |0 > • (9) 

v 2! 



7=0,2,4 

In the same way, one can introduce the n-phonon channel as 

\n >= 4t(4o)10 > • (10) 

If we truncate to the two-phonon states, the corresponding coupling matrix is given by 

/ F(r) \ 

Vcou V = F(r) kuj V2F(r) , (11) 




where, F(r) is defined as /?2/(^)/v47r. 

B. Anharmonic vibrator 

Higher order corrections to Eq. (^) lead to anharmonicities in the vibrational excita- 
tion. The effects of anharmonicity in surface vibrations can be described in many different 
ways. Among them, the interacting boson model (IBM) |50| in the vibrational (U(5)) limit 
provides a convenient calculational framework to discuss the anharmonic effects. The vi- 
brational limit of the IBM and the anharmonic vibrator (AHV) in the geometrical model 
are very similar, the only difference coming from the finite number of bosons in the former 



6], 57]. The eigenvalues of the intrinsic Hamiltonian H int in the U(5) limit are given 

by a 

E ndVl = en d + an d (n d + 4) + (3v(v + 3) + 7/(1 + 1), (12) 

where n d , v, and I are the quantum numbers giving the number of <i-bosons, the <i-boson 
seniority, and the intrinsic angular momentum, respectively, e, a, (3, and 7 are adjustable 
parameters. The first term gives equally spaced and degenerate phonon spectra, while 
the splitting of multi-phonon multiplets due to the anharmonic effects are caused by the 
remaining terms. 

A model for subbarrier fusion reactions, which uses the IBM to describe the effects 
of channel couplings, has been developed by Balantekin et al. |]5~T , |58| , [59[] . They assume 



that the coupling Hamiltonian has a similar form to that of the collective model Eq. (0). 
In the linear coupling and the no-Coriolis approximations, it is given as 

Vmup (r,® = -J=L=f(r)Q 20 . (13) 
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Here, N is the boson number and we have introduced the scaling of the coupling strength 
with y/N to ensure the equivalence of the IBM and the geometric model results in the 
large N limit JSJ, [55]. Q20 is the quadrupole operator in the IBM, which is given by 

Q20 = s j d + sdl + X2{d ] d) i £\ (14) 

where tilde is defined as d^ = (— ) m gL^. Systematic studies of subbarrier fusion reactions 
in this model indicate that the coupling strengths used in this model are very similar to 
those in the geometrical model |]5"8| |. We, therefore, assume that the coupling strength 



in Eq. ([13]) is the same as that in the harmonic limit of the geometrical model, Eq. (g). 
This approximation is valid when the anharmonic effects are not so large, so as to allow 
description of vibrational mode of excitations in terms of the U(5) limit. 

In anharmonic vibrators, phonon states within a given multiplet are split in energy. 
Thus one has to treat each state separately as a different channel. The introduction of the 
multi-phonon channel, as in the harmonic limit, is possible only when the multi-phonon 
states are degenerate in the excitation energies. It will be shown in the next section 
that one can safely neglect the effects of this energy splitting on fusion cross section as 
well as on the fusion barrier distribution. Accordingly, in order to keep the calculations 
simple, we assume in the following that the multi-phonon multiplets are degenerate in 
their excitation energy unless explicitly mentioned. The wave function of the n-phonon 
channel in the framework of the IBM then reads 

\n >= 1 J sY- n (4r\0 >, (15) 

and the corresponding coupling matrix, truncated to the two-phonon states, is given by 



/ F(r) \ 



Vcoup 



F(r) huj - T7 J mX 2F(r) y/2(l -l/N)F(r) 



(16) 



V y/2(l - T/N)F{r) 2hu + 5- ^JfyF{r) ) 



The parameter 5 has been introduced to represent deviation of the energy spectrum 
from the equal spacing in the harmonic oscillator limit. When the X2 parameter in the 
quadrupole operator is zero, quadrupole moments of all states vanish, and one obtains the 
harmonic limit in the large N limit. Non-zero values of \2 generate quadrupole moments 
and, together with finite boson number, are responsible for the anharmonicities in electric 
transitions. The differences in the coupling scheme in the harmonic oscillator limit and 
in the case with anharmonicities are summarized schematically in Fig. 1. 

C. All order couplings 

It has been pointed out that the linear coupling approximation for the nuclear coupling 
is not adequate in heavy- ion fusion reactions [^T[ [52|| . Higher order couplings have to be 
included in order to describe realistic systems and compare the theoretical calculations 
with the experimental data. On the other hand, non-linear terms were shown to be 
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negligible for the Coulomb couplings |52] . Higher order corrections to the nuclear coupling 
were discussed in Ref. |52] in the case of a harmonic oscillator. The coupling matrix 
element between n- and m- phonon channels in the harmonic limit is given by 



nm \ 



dxu* n (x)u r 



-Vn 



1 + exp (r - R - <J^R T x) /a 



(17) 



Here we have adopted the no-Coriolis approximation and the nuclear potential is assumed 
to have a Woods-Saxon form. Note that this equation includes the bare nuclear potential 
Vn in the total Hamiltonian (Eq. ([]])). u n {x) is the wavefunction of the n-th excited state 
of the harmonic oscillator (see Eq. (|T0| ) ) and is given by 



u n {x) 



H n (x/V2a Q ) e~ x2 l^ 



where H n {x) is the Hermite polynomial. 

Extension of the anharmonic oscillator coupling in the linear order approximation to 
all orders can be carried out in a similar manner. Following Refs. [[[1], 58, [59] , we assume 
that the nuclear coupling Hamiltonian is given in the no-Coriolis approximation by 



-v a 



1 + exp (r — R 



Ri 



/47T7V 



Q20) /a 



(19) 



The matrix elements of Eq. (|T^) can be evaluated most easily by introducing the inter- 
action representation which diagonalizes the quadrupole operator Q20 HH- Since phonon 
states \n, I, M) with M ^ do not couple to the ground state in the no-Coriolis ap- 
proximation, we do not need to consider the terms $ m d m with m 7^ in the quadrupole 
operator, Eq. (|Tjj). Accordingly, we diagonalize the operator 



Q20 — s^cIq + sd\ + c^d^do, 
where Co = — ^2/7x2, in the interaction basis as 

Q20 = e + B ] + B + + e_B ] _B_. 



(20) 



(21) 



where B± (B±) are creation (annihilation) operators of the eigenbosons and e± are their 
eigenvalues. After working out the commutators 



Q20, B± 



e±Bl 



we find that they are given by 



4 

Bl 
e± 



y s ] + y 2 dl, 



~V2S 

c ± 



4 + d 



(22) 



(23) 
(24) 

(25) 
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with y and y 2 are denned as 



2/o 



JJ2 



V 4 


+■ Co - c 


2 \ 


A+cg 




f Co + c 


2 \ 





(26) 
(27) 



The eigenvectors of the quadrupole operator are then given by 



A/n + !n_! 



(5- 



10 >, 



:28i 



and the matrix elements of the nuclear coupling Hamiltonian between n- and m- phonon 
channels read 



nm V 



N 

k=0 1 + exp ^ 

k 

i,j=0 







i? — -R2 



(A;e+ + (N - ife)e_)) /a 

'(J)' 



v^ttJV 

N — k 
N — n — i 



N-k 
N — m — j 



x 



yJ(N - n)\n\(N - m)\m\ 
k\(N — k)\ 



2i+2j-2k+m+n 2k+2N-m-n-2i~2j 

y V2 



(29) 



In evaluating Eq. (|29|), we used the fact that the sum of the number of each eigenfunction 
bosons, n + + n_, must be equal to the total boson number N. 

D. Coupling to octupole mode 

In spherical nuclei, the low-lying octupole vibrational state has a large collectivity 
and their excitations play an important role in subbarrier fusion reactions. Therefore, in 
order to apply the models which we discussed above to realistic systems, it is necessary to 
extend them so that they include the octupole mode as well. The coupling Hamiltonian 
in the harmonic limit can be extended in a straightforward way and reads 



-Vn 



a) 



and 



1 + exp (r — i?i 



4 --Rt«20 



l-R T (y 30 ) j a 



(30) 




5 3 0R 3 

a 2 o + -z P Z T e 



2 ±L T_ 

4 \/ 4^ 




"30, 



(31) 



v ' 5 r 6 V 4vr 7 r y 

for the nuclear and the Coulomb couplings, respectively. Here a 2 o and a 30 are the surface 
coordinates for the quadrupole and the octupole vibrations, respectively. In order to 
take into account the octupole mode of excitation in the anharmonic oscillator coupling, 
we use the vibrational limit of the sdf- IBM [ 50" , pH| ]. This model has been successfully 
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used to describe negative parity states in rotational |K| as well as vibrational |6T|, |62], p^j 
nuclei. As a natural extension, we assume the following coupling Hamiltonians based on 
this model for the nuclear and the Coulomb couplings. 

V [N) (r = — (32) 

1 + exp [(r -Ro- Rr^Q« ~ RtJ^Qso) /a] ' 

= ^We^^o.^V^g, (33) 

Here, /?3 is the octupole deformation parameter, and we take the quadrupole and the 
octupole operators in the sdf- IBM as 

Q 2 = Jd+ sS + x 2 (dW 2) + X2/(/ t /) (2) , (34) 
Q 3 = sp + X3W^ + h.c., (35) 

respectively, where / M is defined as (— ) 3+M /- At . We will apply these models in Sec. IV to 
analyze the 16 O + 148 Sm fusion reaction. 

III. ANHARMONIC EFFECTS ON FUSION BARRIER DISTRIBUTIONS 

In order to discuss the effects of anharmonicity on subbarrier fusion reactions, we 
perform in this section a series of model calculations without particular regard for the 
physical values of the parameters. We consider a fictitious nucleus whose excitation 
energies are given by Eq. ( |T2] ) with e = 800, a = 10, (5 = and 7 = 25, all in keV. The 
excitation energy of the first 2 + state of this nucleus is 1 MeV and the double phonon 
states are split by 0.5 MeV. We take the total boson number N to be 4 and thus take into 
account coupling of up to four phonon states. The other parameters used are X2 = _ 3 
(see Eq. fll4|) ), the quadrupole deformation parameter @ 2 — 0.25 and the target radius 
Rt = 6.35 fm. Finally, the Christensen-Winther parameterization of the Woods-Saxon 
potential [Q for the 16 + 148 Sm system is used for the ion-ion potential Vn- 



As we discussed in the previous section, introduction of the multi-phonon channels 
significantly reduces the dimension of the coupled- channels equations. To see whether this 
approximation is good, a calculation is first made by taking fully into account the effects 
of the splittings in the phonon spectrum. In order to keep the calculation simple, we make 
here the linear coupling approximation. In Fig. 2, we show the excitation function of the 
fusion cross section (the upper panel) and the barrier distribution (the lower panel) that 
follow from using the U(5) limit as described above (solid line). Barrier distributions are 
obtained using a point- difference formula |1(J with AE =2 MeV in the laboratory frame. 
To see the effect of the splittings of the multi-phonon states, we have repeated the same 
calculation by assuming that the n-phonon multiplets are degenerate at the excitation 
energy of the first (2n) + state. We found that this prescription leads to no visible change 
in the fusion barrier distribution compared to the actual case with splitt energy levels. In 
the following, we shall therefore assume a degenerate multi-phonon spectrum and use a 
single multi-phonon channel representing all the states in a multiplet. 
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In Fig. 2, we also compare the results for the anharmonic vibrator with those in the 
harmonic limit (the dotted line). One can observe that the anharmonic effects in the 
vibrational excitations lead to a significant change in the barrier distribution, though the 
excitation function of the fusion cross section itself is modified only marginally. 

A. Anharmonicities in excitation energies 

In all the following calculations, we take into account the couplings to all orders. We 
first discuss the effects of the deviation of the energy spectrum from the harmonic limit. 
In order to isolate these effects, we truncate the phonon spectrum at the double-phonon 
states in this subsection. Fig. 3 shows the dependence of the fusion cross section (the 
upper panel) and the fusion barrier distribution (the lower panel) on the parameter 5 in 
Eq. ([16]). The cross section and the barrier distribution are calculated for three different 
values of 5. The solid line is the result when 5 = 0, while the dotted and the dashed 
lines are obtained by setting 5 to be 0.5 and —0.5 MeV, respectively. Despite the fact 
that an unrealistically large value is used for 5 so as to maximize its effect, the barrier 
distribution changes only slightly for different choices of 5. This indicates that the main 
effects of anharmonicity on fusion barrier distributions come from the deviation of the 
transition probabilities from the harmonic limit, including the reorientation effect, rather 
than anharmonicities in the level energies. Note that this observation does not necessarily 
mean that fusion reactions are not sensitive to the excitation energy of the phonon states. 
To demonstrate this, we show in Fig. 3 also the result when both Uuj and 5 is set to be 
zero (the dot-dashed line). One can see sizable effects of the finite excitation energy of the 
phonon states on the fusion barrier distribution as well as fusion cross section. We thus 
conclude that, though the fusion cross section and the barrier distribution are sensitive 
to the energy of the single-phonon state, they do not depend so much on the excitation 
energies of the multi-phonon states once the phonon quanta hu is fixed. 

B. Reorientation effects 



One of the pronounced features of an anharmonic vibrator is that the excited states 



have non-zero quadrupole moments |20], [HJ]. In the IBM, the E2 operator is defined as 
T(E2) = cbQ2, where the quadrupole operator was introduced in Eqs. ( fT^D and fl34|) for 
the sd- and sdf-lBM, respectively, e# being is the effective charge. In the U(5) limit, the 
E2 effective charge es is related to the E2 transition probability by 

B(E2; 0+ -> 2+) = 5e|iV. (36) 

Using the definition of the static quadrupole moment of a state with spin I 



Q(I) = J 1 -^<II\T(E2)\II>, (37) 



we obtain for the quadrupole moment of the 2 + state 



Q{2+) = vlr eBX2 - (3 ^ 



ii 



The corresponding formula for the first 3 state reads 



Q(3~) = \I^BX2f. (39) 

Fig. 4 shows the influence of the sign of the quadrupole moment of the first 2 + state on 
the fusion cross section (the upper panel) and the fusion barrier distribution (the lower 
panel). The solid line corresponds to the negative static quadrupole moment, while the 
dotted line is obtained by inverting the sign of the \2 parameter in Eq. flT4j). As seen 
from Eq. (|3~8D , this is equivalent to taking the opposite sign for the quadrupole moment 
of the excited state. The dashed line is the result when \i — 0. Fig. 4 demonstrates 
that the fusion cross section, and especially fusion barrier distribution, strongly depend 
on the sign of the \2 parameter, hence on that of the static quadrupole moment. This 
fact will be used in the next section to determine the quadrupole moments of the first 2 + 
and 3~ states in 148 Sm from the experimental fusion barrier distribution for the 16 + 
148 Sm reaction. 

C. Finite N effects 

The other important effect of the deviation of the electric transition rates from the 
harmonic limit is caused by the finiteness of the boson number. Eq. ( |i~6D and Fig. 1 indicate 
that the anharmonic effects weaken the coupling between the one- and two-phonon states 
by the factor y 1 — 1/iV compared with that in the harmonic limit. For small values of 
the boson number N, this factor is significantly smaller than one and large anharmonic 
effects on fusion reactions are expected. In order to demonstrate the finite boson number 
effects, we show in Fig. 4 also the results in the harmonic limit by the dot-dashed line. 
One can see a significant difference between this result and that for the zero quadrupole 
moment (the dashed line), indicating the importance of the finite boson number effect. 

Finite boson number effects can also be studied in another, perhaps more instruc- 
tive way. Since the couplings to the multi-phonon states are weaker than those in the 
harmonic limit if the anharmonic effects are present, one expects that fusion excitation 
function converges more rapidly compared with the case in the absence of the anharmonic 
effects. Fig. 5 shows how the fusion cross sections (the upper panel) and the fusion barrier 
distributions (the lower panel) converge with the phonon number for fixed total boson 
number N = A. Although the fusion cross section and the fusion barrier distribution 
obtained by including only the single-phonon excitations are quite different from those 
in the double-phonon couplings, the difference between the two- and the three-phonon 
couplings are small. We observe that the fusion barrier distribution almost converges at 
the three-phonon level. The corresponding calculations in the harmonic limit are per- 
formed for comparison. These results are shown in Fig. 6. Contrary to the results for 
the anharmonic vibrator, the barrier distribution obtained by taking into account up to 
the four-phonon excitations is still significantly different from that obtained by includ- 
ing up to the three-phonon excitations. We thus conclude that the finite boson number 
significantly influences the role of higher phonon states in determining the fusion barrier 
distribution. 
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IV. 16 + 148 Sm REACTION: COMPARISON WITH 
EXPERIMENTAL DATA 



We now apply the formalism described in Sec. II to a realistic system. We analyze in 
particular the 16 O + 148 Sm reaction, whose excitation function of the fusion cross section 
was recently measured with high accuracy [|IIJ. We take both the quadrupole and the 
octupole vibrational excitations in 148 Sm into account, while the excitations of 16 are 
not explicitly included in the coupled-channels calculations. The latter have been shown 
to lead only to the shift of the fusion barrier distribution in energies without significantly 
altering its shape |65j and hence can be incorporated in the choice of the bare potential. 

In order to qualitatively understand the effects of the channel couplings on this reac- 
tion, calculations are first performed by assuming the harmonic limit. The quadrupole and 
the octupole deformation parameters are estimated from the experimental electric tran- 
sition probabilities between the ground and the one-phonon states |66| to be /3 2 = 0.182 
and Ps = 0.236, respectively, by assuming the target radius of Rt = 1.06 A 1 / 3 fm. The 
excitation energies of the phonon states in 148 Sm are 0.55 and 1.16 MeV, for the first 2 + 
and 3~ states, respectively. The parameters for the bare ion-ion potential are obtained 
by fitting the experimental fusion cross sections. It has been pointed out that the ef- 
fects of channel couplings play an important role in determining the bare potential for 
the 16 + 144 Sm reaction p7| . Accordingly, we fit the experimental data by assuming 
the three phonon couplings in the harmonic limit (see the following discussion). We use 
the experimental data between 200 mb and 400 mb to determine the bare potential. We 
choose this range because, more details of the channel couplings would be important in 
the lower energy region, while some other effects, e.g. the angular momentum truncation 
or the dissipation mechanism might play some role [^] in the higher energy region. The 
best fit parameters which we obtain are Vq = 155.1 MeV, Rq = 0.95(Ap 3 + A^ 3 ) fm, and 
a = 1.05 fm, for the depth, the radius, and the diffuseness parameters of the Woods-Saxon 
potential, respectively. Note that the experimental data at high energy region still require 
a large surface diffuseness parameter fl(| even after including the effects of the channel 
couplings. The origin of this is not fully understood. 

The results for the fusion barrier distribution for the 16 O + 148 Sm reaction are shown 
in Fig. 7. The panels differ from each other by the number of the quadrupole phonons 



coupled. The experimental data are taken from ref. ||10|| . The point difference formula 



with AE = 2 MeV in the laboratory frame is used to obtain the fusion barrier distribution. 
Calculations, which assume that the octupole vibration in 148 Sm has only a single-phonon 
state, are indicated by the dotted lines in Fig. 7. In the three phonon calculation, for ex- 
ample, we included the 2+, 3" , 2+ (g) 2+ 2+ (g) 3~, 2+ (g) 2+ (g) 2+ and 2+ (g) 2+ (g> 3~ states 
in 148 Sm in the harmonic limit. The coupling scheme in the other panels is defined in the 
same way. All the dotted lines which are obtained by assuming the single-octupole-phonon 
excitations in 148 Sm do not account for the experimental fusion barrier distribution. How- 
ever, double-phonon excitations have been found in the neighbouring nuclei, i.e. 144 Sm 
39 , 40~1) |4l|1 , 146 Sm ||43|| , and 147 Sm [42]|. We therefore repeat the same calculations by 



assuming the double-octupole excitations in 148 Sm (the solid lines). One can now observe 
that the experimental data can be explained reasonably well by the three-phonon calcu- 
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lations. The fit is not as good if the four phonon quadrupole excitations are included in 
the coupled channels analysis. These calculations in the harmonic limit thus suggest that 
there are strong couplings up to the three quadrupole phonon states in 148 Sm and the 
coupling between the three and four phonon states are not as strong as that expected from 
the harmonic oscillator coupling. They also suggest that there might be double-octupole 
phonon excitations in 148 Sm. 

Calculations which take into account the anharmonicities of the vibrational excitations 
are performed next. Following the results of the analysis in the harmonic oscillator limit, 
we include the octupole excitations up to the double phonon levels and take the total 
boson number N = 4. The latter is consistent with the existence of the Z = 64 subshell 



closure |56j]. The standard prescription for the boson number (i.e. counting pairs of 
nucleons above or below the nearest shell closure) would give N = 8 for 148 Sm. However, 
it is well known that due to the Z = 64 subshell closure, the effective boson number is 
much smaller. The suggested effective number in the literature for the proton bosons 



varies between N w = 1 and 3 |69], fTOH , leading to the total boson number to be between 
3 and 5 for 148 Sm. There are three other parameters, i.e. \2i X2f, an d Xs parameters in 
the transition operators Eqs. fl3~4|) and fl35|), which need to be determined. Two of them, 
X2 and X2f, are directly related to the quadrupole moment of the phonon states, and 
the other the coupling between the quadrupole and the octupole modes. Though there 



exist experimental data on the quadrupole moment of the first 2 + state of 148 Sm [ 66H , the 
experimental uncertainty is still large. There are no other experimental data which can be 
used to fix the \2f and \z parameters. Therefore, a \ 2 fit to the experimental fusion cross 
sections is carried to determine all the three parameters. In the fitting procedure, all the 
experimental data below 200 mb are used except for the one at E; a b= 69.36 MeV, which 
lies outside the systematics. The best fit values are X2 = —3.12 ± 0.77, X2f = 4.63 ± 0.43, 
and X3 — —1-99 ± 0.26. The solid line in Fig. 8 shows the fusion cross section (the upper 
panel) and the fusion barrier distribution (the lower panel) thus obtained. Following the 
discussion in the previous section, harmonic spectra for the excitation energies of the 
phonon states are assumed. One observes that the shape of the fusion barrier distribution 
obtained by including up to four phonon states with anharmonicities is very similar to 
that for the three phonon couplings in the harmonic limit (see Fig. 7) and that they agree 
quite well with the experimental data. This situation is similar to that observed in ref. 
|46|1 , where it has been shown that the shape of the fusion barrier distribution for the 16 O 
+ 144 Sm reaction obtained by including up to single phonon states in the harmonic limit 
is very similar to that for the double phonon couplings with anharmonicities. 

In the last section, we discussed two main effects of the anharmonicities in the electric 
transitions, i.e. the reorientation and the finite N effects on subbarrier fusion reactions. In 
order to see each effect separately, we repeat the same calculation by setting both the \2 
and X2f parameters to be zero. As was discussed in the previous section, this prescription 
is equivalent to taking only the finite N effects into account. The result is shown by 
the dashed line in Fig. 8. The figure also shows for comparison the result obtained by 
assuming that there are four quadrupole phonon states in the harmonic oscillator limit 
in 148 Sm (the dotted line). The difference between the dotted and the dashed lines, and 
also between the dashed and the solid lines, is significant. This suggests that both the 
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reorientation and the finite N effects play an important role in the fusion reaction. 

From the \2 and \2f parameters which are obtained by the x 2 fit to the experimental 
fusion cross sections, we estimate the quadrupole moments of the first 2 + and 3 _ states 
of 148 Sm. The E2 effective charge e# is estimated to be 0.19 eb from the experimental 
B(E2) value (see Eq. @). Using Eqs. ([38]) and @, they are estimated to be Q(2+) = 
—1.00 ± 0.25 b and Q(3~) = +1.52 ± 0.14 b, respectively. The quadrupole moment of the 
first 2 + state thus obtained is very close to that measured from the Coulomb excitation 



technique, i.e. —0.97 ± 0.276 [66]. There is no experimental data to compare our result 
for the quadrupole moment of the first 3~ state. We can, however, test the consistency of 
the fit by taking its opposite sign in the coupled-channels calculations. Fig. 9 shows the 
sensitivity of the fusion cross section and the fusion barrier distribution to the sign of the 
quadrupole moment of the first 3~ state. The solid line corresponds to the optimal choice 
for the sign of the first 3~ state, while the dotted line is obtained by inverting it. The use 
of the incorrect sign of the quadrupole moment destroys the good fit to the experimental 
data. Our calculations thus strongly suggest that the heavy-ion fusion reactions at low 
energies can provide an alternative method to determine the sign as well as the magnitude 
of the quadrupole moments of the low-lying excited states in spherical nuclei. 

V. SUMMARY AND CONCLUDING REMARKS 

We discussed the effects of anharmonic phonon excitations on heavy-ion fusion reac- 
tions at subbarrier energies and showed that they play an important role. We showed 
that the vibrational limit of the interacting boson model provides a useful framework to 
address these questions. There are mainly three effects of anharmonicities; the anhar- 
monicity in the excitation energy spectra, the reorientation effects, and the finite boson 
number effects in the transition strength. We showed that the anharmonic effects in the 
energy spectra play only a minor role and the main effects come from the deviation of 
transition probabilities from the harmonic limit, i.e. both the reorientation and the finite 
N effects. We found that the fusion barrier distribution strongly depends on the sign 
of the quadrupole moment of excited states. Using these properties, we have analyzed 
the 16 O + 148 Sm fusion reaction and discussed the anharmonic properties of the phonon 
excitations in 148 Sm. It was found that the best fit to the experimental data requires 
a negative and a positive quadrupole moments for the first 2 + and 3~ states of 148 Sm, 
respectively. Since the quadrupole moment of the first 3~ state of 144 Sm has been found 



to be negative |[46|| , it would be interesting to measure that of the nucleus in between, i.e. 
146 Sm. 

We have shown that calculations in the harmonic limit provide only qualitative re- 
sults and the realistic situations are much more complex due to anharmonicities. Since 
harmonic calculations suggested that a fusion barrier distribution is quite sensitive to the 
number of phonon states excited during the fusion, it has been hoped that subbarrier 
fusion reactions offer an alternative method to identify the existence of higher phonon 
states. Actually, some analyses, aiming at identifying multi-phonon states, have been 
carried out. It is, however, apparent from our results that care has to be taken in such 
analyses. Our calculations show that fusion barrier distribution converges more rapidly 
with the number of phonons when the anharmonic effects are present compared with the 



15 



case in their absence. Harmonic calculations may thus underestimate the maximum num- 
ber of phonons in vibrational nuclei. It is vital to take into account anharmonic effects in 
order to properly identify the maximum number of phonons. 
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Figure Captions 



Fig. 1: Schematic representation of the coupling scheme in vibrational models. The up- 
per part is for the anharmonic vibrator, while the lower part for the harmonic oscillator. 
See text for the notation. 



Fig. 2: Effects of anharmonicity of the quadrupole vibration on the excitation function 
of the fusion cross section (the upper panel) and fusion barrier distribution (the lower 
panel). The dotted line is the result in the harmonic limit, while the solid line takes the 
anharmonic effects of the vibrational excitations into account. The linear order coupling 
is assumed. 



Fig. 3: Dependence of the fusion cross section (the upper panel) and fusion barrier 
distribution (the lower panel) on the deviation of the energy spectrum from the harmonic 
limit. The solid, the dotted, and the dashed lines are obtained by setting the parameter S 
in Eq. flTBD to be zero, 0.5, and —0.5 MeV, respectively. The dot-dashed line is the result 
in the degenerate limit in the anharmonic vibrator coupling. The phonon spectrum is 
truncated at the double phonon levels. The full order couplings are included. 

Fig. 4: Dependence of the fusion cross section (the upper panel) and fusion barrier 
distribution (the lower panel) on the sign of the quadrupole moment of the first 2 + state. 
The solid and the dotted lines correspond to the cases for the negative and the positive 
static quadrupole moments, respectively. The dashed line was obtained by setting the 
quadrupole moment of the excited states to be zero. The dot-dashed line is the result of 
the corresponding calculations in the harmonic limit. 

Fig. 5: Convergence of the fusion cross section (the upper panel) and fusion barrier 
distribution (the lower panel) as functions of the number of phonon states included in the 
coupled-channels equations, which are indicated in the inset. The total boson number N 
is fixed to be four in all cases. 



Fig. 6: The same as Fig. 5, but for the harmonic oscillator coupling. 



Fig. 7: Comparison of the experimental fusion barrier distribution with the results of 
coupled-channels calculations in the harmonic limit for the 16 + 148 Sm reaction. The 
dotted lines include the single octupole excitations, while the solid lines take the double 
octupole phonon couplings into account. The experimental data are taken from Ref. (TO . 



Fig. 8: Comparison of the experimental fusion cross section (the upper panel) and fusion 
barrier distribution (the lower panel) with the coupled-channels calculations for the 16 
_l_ i48g m reac ti on The experimental data are taken from Ref. The anharmonic 

effects in the quadrupole and the octupole vibrational excitations in 148 Sm are taken into 
account in the sdf- IBM. The solid line is the result when the best fit parameters are used, 
while the dashed line is obtained by setting the quadrupole moments of all the excited 
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states to be zero. The results including the four phonon couplings in the harmonic limit 
is denoted by the dotted line. 

Fig. 9: Dependence of the fusion cross section (the upper panel) and fusion barrier 
distribution (the lower panel) on the sign of the static quadrupole moment of the first 
3~ state for the 16 + 148 Sm reaction. The solid and the dotted lines correspond to the 
cases where the quadrupole moment of the first 3— state in 148 Sm is positive and negative, 
respectively. 
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